#ifndef _MATRIXMATH_
#define _MATRIXMATH_

#include <vl/VLd.h>
#include <vector>
using namespace std;

namespace moremath{
	Matd getInnerProdMat(vector<Vecd> _a, vector<Vecd> _b);
	Matd getInnerProdMat(Matd _a, Matd _b);
	void getColMeanCenteredMat(Matd &_m, Matd *_newm, Vecd *_mean);
	double traceMult(Matd &_A, Matd &_B);	// trace of A*B
	double traceMultT(Matd &_A, Matd &_B);	// trace of A*B^T
}

#endif
